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Abstract 

We consider anew some puzzling aspects of the equivalence of the quantum field theoretical description of 
Bremsstrahlung from the inertial and accelerated observer's perspectives. More concretely, we focus on the 
seemingly paradoxical situation that arises when noting that the radiating source is in thermal equilibrium 
with the thermal state of the quantum field in the wedge in which it is located, and thus its presence does 
not change there the state of the field, while it clearly does not affect the state of the field on the opposite 
wedge. How then is the state of the quantum field on the future wedge changed, as it must in order to 
account for the changed energy momentum tensor there? This and related issues are carefully discussed. 
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I. INTRODUCTION 



The topic of radiation by uniformly accelerated charges has often been the source of much 
puzzlement and confusion, particularly when considered in the light of the equivalence principle. 
Much of the confusion has been removed by the realization that for observers coaccelerating with 
the charge there are regions of spacetime that are inaccessible. In effect, in the classical context, it 
has been shown that for an electromagnetic charge in uniform acceleration, the classical radiation 
field, as described by accelerating observers, is zero at every point in the region that is accessible 
to them (known as the Rindler wedge, R, see Fig. I) [1]. This would then remove any apparent 
contradiction between the EP and the known Bremsstrahlung. For both the inertial and accelerated 
observers there is radiation but, for the accelerated ones, such radiation lies beyond the regime 
where the static description is valid. 

In the quantum version of this situation the question is posed in terms of emission of pho- 
tons rather than the evaluation of radiation fields. In fact, it was shown that the standard 
Bremsstrahlung when viewed from the point of view of the accelerated observers -a point of view 
called Rindler quantization- acquires a very particular interpretation. Actually, as will be ex- 
plained in more detail below, the coincidence for the prediction of photon emission rates between 
the inertial and accelerated descriptions makes fundamental use of the Unruh effect [2], which 
states that from an accelerated frame comoving with the charge, the standard Minkowskian vac- 
uum state corresponds to a thermal state. Moreover, it is well known that for a detector uniformly 
accelerating in the inertial vacuum the process for which the detector absorbs a particle from the 
bath (a Rindler particle) as seen by a comoving (accelerating) observer is equivalent to the emission 
from the detector of a Minkowski particle as seen by an inertial observer [3] . Then from the point 
of view of an inertial observer the accelerating charge emits particles while, from an accelerated 
viewpoint, the charge -which is static- will emit and absorb Rindler particles to and from the 
bath. 

The restriction of this effect to wedge R has been analyzed in Ref. [4], where it is shown that the 
emission rate of photons with fixed transverse momentum in the inertial frame coincides with the 
combined rate of emission and absorption of zero-energy Rindler photons with the same transverse 
momentum in the accelerated frame. Thus, this result gives a clear notion of the the physical 
equivalence between inertial and accelerated descriptions of Bremsstrahlung. However, and as is 
often the case in this field, the answer to one question brings in further puzzlement, and the need 
to answer a further one. 

The calculation mentioned above, makes fundamental use of the so called zero energy Rindler 
modes. The reason is that for accelerated observers the charge is static and thus it can only couple 
to modes of zero frequency with respect to Rindler time. In fact, due to the expression of the form 
x oo that appears in calculations involving zero energy particles in the analysis of Ref. [4], it is 
required the introduction, for the purpose of regularization, of a small frequency of oscillation 6 
for the source, which allows one to work with finite energy modes and which at the end is taken to 
0. In that work the authors considered also the question of whether or not an accelerated observer 
sees any difference in the thermal bath due to the emission and absorption of Rindler photons. 
They note that in the limit of zero energy the transition rate from an n-photon state to an (n + 1)- 
photon state and the rate of the inverse process became equal, thus the accelerated charge leaves 
the thermal bath undisrupted. From this one can conclude that the source is in thermal equilibrium 
with the quantum field. Then, from the point of view of an accelerated observer, in R, there will 
be no difference between the initial state of the field (the initial thermal state) and the state 
generated by the interaction with the accelerating charge (note the analogy with the description 
of the situation in the classical context that we explained above) . Similarly, the state of the field 
in the second Rindler wedge will remain the initial thermal state, as that region could not possibly 
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be, in anyway, influenced by what is going on in a causally disconnected region of spacetime. We 
could think that this conflicts with the fact that EPR influences are allowed but this is misleading: 
when talking about the state in the L Rindler wedge we mean the corresponding density matrix 
and this could not change due to actions taken on the right wedge for otherwise there would be 
operators pertaining to the left wedge whose expectation values could be used to determine what 
has occured in the right Rindler wedge. Nevertheless, it is clear that in wedge F (see Fig. I) 
there would be a detectable change in the state of the field and, in particular, the expectation of 
the energy momentum tensor in this region should be different from what it would have been if 
there was no interaction with the charge. In effect, this change could be computed in an inertial 
quantization of the field and would correspond to the final state containing the Minkowski photons 
emitted by the the accelerated charge. On the other hand, the quantum description of the state of 
the field as seen by the accelerated observer is given by the Unruh quantization scheme [2, 5]. In 
this description, the restriction to the L and R wedges of both, the Minkowski vacuum state and 
the state resulting from the interaction with the charge is, in both thermal bath, and thus 

it seems neither could contain the information regarding what has changed in F. The issue is then 
whether this situation can be analyzed in the language appropriate to accelerated observers and 
how would, in that case, the information about the changed situation in F be codified? 

The fact that Unruh modes can be expressed by superpositions of specific modes that extend 
distributionally to the whole spacetime (known as boost modes, see Refs. [6, 7]), opens the door to 
analysis of issues related to physical questions outside the double wedge, and in particular in wedge 
F, carried out in the language appropriate for accelerated observers. In this work, we embark on 
such analysis for the case of a scalar field and an accelerated scalar source, and analyze from the 
new perspective the change in the expectation value of the energy-momentum tensor operator when 
evaluated at points in wedge F. 

States in the Unruh quantization can be seen as states of a composite system where each of 
the components is the quantum field restricted to either wedges L or R. Then, having the density 
matrix p for a state, one can describe the physics, for example, in wedge L (R) by tracing out 
in p the right (left) degrees of freedom. From this procedure we obtain a density matrix p L (p R ) 
describing a state in wedge L (R). It is well known that when a field observable A is localized in 
either L or R, then the expectation value of A is determined completely by p L or p R respectively. 
On the other hand, when A is localized out of the double wedge, it is clear that its expectation 
value would not, in general, be determined solely by information encoded in p L or p R . 

There should exist some object containing this extra information carried by the state which 
controls how do left and right parts combine which is the extra element necessary for describing 
completely the state in all of Minkowski spacetime. We identify this object as the entanglement 
matrix. In the case of interest we were concerned about the expectation value of the energy 
momentum tensor in the future wedge and it seemed natural to expect its change to be encoded 
in the change of the entanglement matrix. This was our initial assumption, and the work intended 
to see how exactly is such information encoded in this case. As we will see this expectation was 
mistaken and the answer in our case lies, surprisingly, in p R . Precisely how does this matrix codifies 
this information will be elucidated through the rest of the manuscript. 

Regarding the foundational basis of this work we must point out that the issue of extending the 
Unruh quantization outside the double wedge is somewhat subtle. What is very well established in 
the literature is the restriction of global states to the double wedge -for example, the restriction of 
states to the wedge R is interpreted as the state seen by accelerated observers which have available 
only the wedge R— but, as far as we know, the extension of a state in the Unruh quantization to the 
whole MS does not seem to be, at this time, fully investigated. One important point concerning this 
question is that Unruh modes are highly singular at the asymptotes which are taken to coincide 
with the horizon of the particular construction. Thus, the initial data from which one builds up 
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the (one particle) Hilbert space in the Unruh quantization has restrictions at these asymptotes 
which can influence the solutions of the field equation in wedges F and P (this possibility has been 
pointed out in Ref. [7], [6]). As we explained above, in this paper we extend the usage of the Unruh 
description of the field to address questions related to the physics in the wedge F. We do this in 
the way that we consider is the most natural and we find results that are physically consistent and 
which coincide with what one would obtain working in the standard Minkowski quantization of the 
field as we show explicitly in Appendix C. 

The paper is organized as follows. In Sec. II we review the main ideas of the Unruh quantization 
and specify our notation. In Sec. Ill we explain our strategy for analyzing the encoding of the 
information in the final state of the field. We propose a particular decomposition of the density 
operator of the state and use it to write an expression for the change in the expectation value of 
T M „. Nevertheless, the particular form of the source representing an accelerating particle is needed 
in order to obtain an explicit expression for the final state and its density operator. In Sec. IV we 
work out these ideas introducing a scalar accelerating pointlike source with a particular regulator 
and build up the S-matrix operator. In Sect. V we make the final calculations in order to obtain 
the change in {T^ v ) and we evaluate it in wedges R and F. From these results one can see how 
the information of the physical change in F is encoded. In Sec. VI we analyze the case of two 
different sources accelerating in wedges L and R in order to get more insight of the behaviour 
of the density operator of the respective final state. Finally, we end with some discussions and 
the interpretation of the results in Sec. VII. In order to carry out these steps it is important to 
have explicit expressions of the spacetime behaviour of the Unruh modes, which we present in 
Appendix A. The calculations leading to our final results are, though straightforward, very long, 
in Appendix B we give an explicit derivation of the main formula in this work. In Appendix C we 
make the same calculations as in Sec. V for the change in wedge F but using, instead, the standard 
plane wave representation of the field and show that both results coincide. 

To eliminate unnecessary notation we shall work in two dimensional Minkowski space, this will 
allow us to be clearer without losing physical insight. As a matter of fact, working with a scalar 
field in 2D with m / is operationally equivalent to work with the same field in 4D and fixed 
traverse momentum k\ = k 2 + k 2 , with the identification m 2 — > k\ + m 2 . In all this work we shall 
work in units in which c = h = 1 and signature — h 



II. UNRUH QUANTIZATION 

All the comoving observers to a particle moving with uniform proper acceleration a = (a^a 11 ) 1 ^ 2 
have world lines of the form 

t = Csinh(ar) z = Ccosh(ar), (H-l) 

where < ( < oo and r is the proper time of the observer, — oo < r < oo. Eq. (II. 1) can be used 
to give coordinates r, ( to wedge R, which is known as Rindler spacetime. In these coordinates, 
the Minkowski metric becomes 

ds 2 = -Q 2 dT 2 + dQ 2 . {11.2) 

As can be seen from Eq.(II.2), wedge R is a static, globally hyperbolic spacetime so one can use 
standard methods of frequency splitting [5, 8, 9] for quantizing the field in that region. As a first 
step we shall construct directly the one particle Hilbert space for a Klein Gordon field with mass 
m in wedge R, using as timelike Killing field r a ; this is called the Fulling- Rindler quantization [10]. 
Consider the massive Klein Gordon (KG) field equation 

(□-m 2 )0 = O (II.3) 
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restricted to wedge R and take solutions which are positive frequency w.r.t. the Rindler time r 
(in this case (d/dr)^ = b 11 is the time translation generator, cf. Eq. (A.l)) and which vanish 
asymptotically. These are superpositions of the following set of modes 



MO = e-^K^mC), u> > 0, (II.4) 

7T 

where £ = (r, £), and K^x) is the modified Bessel function of the third kind or Macdonald 
function. The normalization factor is chosen such that these modes are 5 normalized w.r.t. the KG 
product: 

(MM kg = \\ tyZVpM - M^M^ = s(u - J) (H.5) 

where is a Cauchy surface for Rindler spacetime. Functions K^imC,) have an essential singu- 
larity (it oscillates "infinitely") in the limit Q — ► [11], so the ipu(£) are not defined at the horizon 
(t — > ±oo and £ — > 0) when w / 0. Modes given by Eq. (II. 4), called Fulling modes, and their 
conjugates ?/>,*,(£) form a complete set of the space of solutions to Eq.(II.3) in Rindler spacetime. 
One constructs the one particle Hilbert space J%r for a quantization of the KG field in R by Cauchy 
completing the space spanned by the ipu{Q (positive frequencies) in the inner product given by 
Eq. (II. 5) (a rigorous construction of this space is given in [12]). The space of states of the field for 
this quantization is the Fock space of J%r, ^(J%r). The field operator in this quantization takes 
the form 

MO = J" (M0?u + MOrt) dw, (II.6) 




Figure 1: The geometry of Rindler space. 
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where the creation and annihilation operators, and r w , satisfy canonical commutation relations. 
Repeating the Fulling-Rindler quantization procedure in wedge L with the L-time translation 
generator given by — b^ (see Eq. (A.l)) and defining the KG product on a Cauchy surface of 
L, one obtains the one particle Hilbert space for the KG field in L, J%l, and a left field operator 
analogous to Eq. (II. 6). 

Wedges L and R, are causally disconnected and thus neither JF(J%l) nor ^■(M' R ) can be used 
to represent the field algebra in Minkowski spacetime, as the Fock space of the standard plane 
wave quantization does. However, one can consider a second quantum field construction in all of 
Minkowski spacetime with one particle Hilbert space given by 

3%j = je> L ® M'r, {11.7) 

and the space of states given by [5] 

jr(j^) = &(jtr R ) ® &{j? l ) . (ii.8) 

The quantum field construction defined by Eq. (II. 7) is called Unruh quantization. The field 
operator in this quantum construction takes the form 

4>{x) = 4> L (x)®i R + i L ®4> R (x). (n.9) 

Whenever there is no confussion, we denote by 4>l(x) the operator 4>l(x) <8> lj? and the analogous 
for 4> R . Clearly, these operators commute, 

[4> L (x),Mx')]=o, (ii.io) 

reflecting the fact that the regions L and R are causally disconnected. 

Constructed in this way, it may seem that in the Unruh quantization the field operator is made 
up of modes which have support only in the double wedge L U R and thus that cannot describe 
the physics outside this region. Nevertheless, recall that the one particle Hilbert spaces J%l and 
J4?r are made up of positive frequency modes which have initial data on either Sl or T, R . These 
two latter Cauchy surfaces can be seen as the restriction of some Minkowski spacetime Cauchy 
surface S to L or R respectively, then the initial data of modes in the left or right quantizations 
of the field define unique solutions of the KG equation in all of Minkowski spacetime. In fact, as 
we explain in the next subsection and in Appendix A, there exist superpositions of plane waves, 
called Unruh modes, which coincide with Fulling modes when restricted to either wedge L or R 
and that are zero when restricted to the opposite wedge [6]. The elements of 3^\j describe modes 
of the field defined distributionally in the whole of Minkowski spacetime. See the last paragraph 
of the next subsection for further discussion. 

A. Boost modes and Unruh modes 

The original quantization approach used by Unruh [2] does not give explicitly the functional 
form of Unruh modes since he works only with their restrictions to the horizons. However, for the 
purposes of this work we do need the functional form in all of MS of these set of modes, which are 
linear combinations of boost modes. We will introduce both sets of modes in this section following 
Ref. [6]. 

Consider the space of global classical solutions to Eq. (II. 3). One can look for solutions in 
this space which have positive frequency w.r.t. the boost parameter r in wedges L and R. This 
class of solutions, Minkowski Bessel Modes were originally presented and studied by Gerlach [13]. 
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Narozhny et.al. [6] and Fulling and Unruh [7] call them boost modes (BM), name that we also 
prefer. Take an (unnormalized) plane wave solution to Eq. (II. 3) of the form 

If{x) = (2tt)- 1 /2 e TiK*- P2 ) ) (n.ll) 

where x = (t, z), uj p = \/p 2 + m 2 > and the upper (lower) sign corresponds to positive (negative) 
frequency w.r.t. inertial time. Changing the variable p to the rapidity 8: 

msinh(6)=p, m cosh(#) = u p , — oo < 6 < oo, (H-12) 

then one can define boost modes as the following superposition of plane waves [6] 

J poo 



, poo 

/ =Fim(t cosh d—z sinh 9) —iuiO in 



(11.13) 



where — oo < to < oo. 

Note that in general, Eq. (11.13) cannot be interpreted as the definition of a function. In 
particular, it is divergent at the origin, (t = 0, z = 0) and thus, cannot stand as a global solution 
to the KG equation if it is considered as a function. However, one can avoid these problems 
if one considers these quantities as distributions, thus requiring the smearing with suitable test 
functions. Therefore, we will consider this set of modes, as well as Unruh modes (cf. Eqs. (11.18)), 
as distributions when constructing the quantum field. Note that this is a consistent procedure 
as the field <fi(x) has, by itself, a distributional character. For the purposes of this work, we will 
consider only test functions of compact support in M. 

Formally, boost modes are orthogonal in the KG inner product (Eq. (II. 5) over Em), 

(B±\B±) KG = ±8(u-n), (11.14) 
and they are eigenfunctions of the boost operator [13] (inside the four wedges of MS), 

O r\ o 

BS± = -iwB± where B = t—+z- = —, (11.15) 

and thus, the modes B^j are positive (negative) frequency KG solutions w.r.t. the boost parameter 
r in wedges L and R whenever to > (lo < 0). Note that in dividing the modes into positive 
and negative frequency modes one is dropping out the u> = mode which could, in principle, be a 
source of trouble [6]. 

From now on we will only use boost modes which are positive frequency w.r.t. inertial time t, 
.£?+, and will drop the superscript: B w = B+. One can use this set of modes as a basis of the 
positive inertial frequency solution space of eq. (II. 3). Actually, one can give a quantization of the 
KG field unitary equivalent to the standard positive frequency plane waves one [6] : 



[x 



/oo 
dw (B„(x)bu + B*(x)bi), (11.16) 
-oo 



where [6 W , b^,] = 5{uj — uj'). Because the transformation is unitary these quantum field descriptions 
share the same vacuum 

k |0m) = a p \0 M ) = 0. (11.17) 
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Unruh's idea [2] for giving a field quantization associated to accelerated observers consists in 
constructing a set of modes R u and L u from combinations of boost modes and their conjugates 
such that they are positive frequency w.r.t. the boost parameter (accelerated time) in wedges R 
and L respectively and zero on the opposite wedge. These modes are defined by 

RUx) = . * . \e™' 2 B u (x) - e-^BlJx)} , (II.18a) 
Y / 2smh(7ro;) L J 

L u {x) = y=J= -. [e^XW " e-™l 2 Bl{x)] , (II.18b) 
y2sinh(7ro;) L J 

where u > 0. Note that, since they are defined in terms of boost modes, these definitions are 
global. These equations can be inverted: 

B u (x) = - L= -. [e™' 2 R„{x) + e-™l 2 Ll{x)] , (II.19a) 
a/2 smh(7ro;) 

B^{x) = 1 [e™/ 2 LUx) + e~™/ 2 R*Jx)] . (II.19b) 

\J I sinn(7ra;J 

It can be seen that, when restricted to R, Unruh modes R w coincide with Fulling modes, 
Eq. (II. 4), and the analogous situation for L w in L [13] (see Appendix A) and, therefore, initial 
data of Unruh modes coincides with that of Fulling modes. It is known that solutions of the field 
equation can be represented by their initial data [8] and thus Cauchy completing the space spanned 
by the L w and R& modes one would obtain, respectively, J^r and ,3%l. 

The field operator in the Unruh quantization takes the form of Eq. (II. 9) with 4>l(x), 4>r(x) 
expressed in terms of Unruh modes: 

/•CO 

<j>(x) = <j> L {x) + $r{x) = I dw(R u ,(x)l L ®r ul + L w (z)f w <g> i fl + H.C.), (11.20) 

Jo 

where f w and l w are annihilation operators in J?(J£Jj) and ^{M'i) respectively. 

We note that recently there has been some controversy about the Unruh quantization. Narozhny 
et.al. [6] claim that the Unruh quantization is not a valid quantization scheme for all of MS (as 
the boost modes quantization is, see Eq. (11.16)). To support this they argue, in particular, that 
the expansion of <f> in Unruh modes, Eq. (11.20), does not exhaust all the degrees of freedom of the 
quantum field in Minkowski space. This is so, they say, because when "evaluating at the origin", 
boost modes have a singularity when to = 0. The authors of Ref. [6] logic is that, consequently, 
in the transition from Eq. (11.16) to Eq. (11.20) — using Eqs. (11.19) — one should take the Cauchy 
principal value of the integral at the origin, which excludes the to = mode. They assert that 
without this mode, the remaining set of boost modes loses the property of being complete and 
then lacks the possibility of spanning every (one particle) state of the field. (Ref. [6], p. 025004-12) 
On the other hand, Fulling and Unruh [7] have argued against this claim. They state that since 
the mode expansion is an integral (Lebesgue measure) and thus one mode is of zero measure, the 
omission of the u = mode is quite harmless in the mode expansion of the field and the Unruh 
quantization is valid for expressing the restrictions to L U R of global MS states. In this work 
we adhere to the position of Fulling and Unruh for the following reasons. As we have said above, 
boost modes should be considered as distributions and thus, evaluating them at one single point 
has no meaning. Moreover, we recall that the field itself evaluated at one point has no meaning 
either, the quantum field is a distributional object and only its convolution with a test function is 
defined. 

In their reply, Fulling and Unruh forcefully argue that the Unruh quantization scheme is valid 
on the double wedge L U R but, however, they indicated that they are not fully confident on the 
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possibility of extending an Unruh state to all of MS ([7] p. 04870 1-2). Let's be more precise on 
this issue. The initial data for the Unruh quantization consists of smooth functions of compact 
support on both Cauchy surfaces and El. For definiteness, take = {(0, z)\z > 0} and 
= {(0,^)1^ < 0}. For these Cauchy surfaces, functions of compact support would be zero at 
the origin (0, 0). Nevertheless, a full MS quantization should consider the set of all Cq° initial data 
over the Cauchy surface 

S M = S L US R U{0}, (11.21) 

which includes, of course, functions which are not zero at the origin. In this respect, Fulling and 
Unruh [7] note that "The treatment of initial data at the origin is mathematically subtle, and data 
at that point may influence the solution of the field equation in regions F and P" . This issue is 
quite relevant for our work. Although we are not giving any formal proof of the fact that the Unruh 
quantization can be extended to all of MS, we show that at least for the particular case we study, 
the Unruh quantization provides the same physical results as the standard flat space quantization. 
This we take as an indication that both quantum descriptions are generally equivalent. 

B. Minkowski vacuum, general states and L-R correlations 

To find a basis of &(3tfc), Eq. (II.8), we need a basis for &(Jtr R ) and &(Jf? L ). For &(Jtr R ) 
we choose an orthonormal basis whose elements are instates \J) R which have a definite number 
of Rindler particles nj. Let J Wm be the number of particles in this state whose frequencies are 
centered in the particular mode uj m , m = 0, 1, . . . . Then, the state | J) R is defined by the set 

oo 

J = { Ju i JbJl ■>•••■> JlOm 1 ■ ■ ■ } ) ^ ] = n J 1 (11.22) 

m=0 

Note that only a finite number of the J Um are / 0. State \ J)r can be build up from the Rindler 
vacuum in the following manner 

\J)R = |^o J U1 ■■■)r = Nj (ftj^o . . . (ftj-W • • • |0) fl , (11.23) 

where Nj is a normalization factor such that (J\J) = 1, and f\ 0m creates a Rindler particle with 
frequency centered at oj m . State \J)r has Rindler energy 

oo 

£(J) = J> m J Wm . (11.24) 

m=0 

This is the state's energy associated with the boost Killing field. Given two different elements \J)r, 
\K) r of this basis we have that 

R (J\K) R = 8{J,K) = 8j ua , Kua ---5j Um jc Um --- , (11.25) 

where 8j^ m ,K^ m are Kronecker deltas. One can obtain a basis of ^■(M'l) in an analogous way. The 
set of all \ J) L ® \K) R where J and K are of the form of Eq. (11.22) is a basis of ^(Jtfjj). Any 
state \g) € ^(Mjj) can be casted in terms of this basis, and it is defined by a particular function 
G(J,K): 

\g) = Y J G(J,K)\J) L \K) R . (11.26) 

J,K 
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Here, the sums run over the space of all possible particle distributions J and K, that is, they are 
of the form 

oo oo 
J m=0 K ldm =Q 

Entangled states in the Unruh quantization (those which cannot be expressed in the form \g\)L <8> 
\g2)R) have non trivial correlations between left and right states. The inertial Minkowski vacuum, 
\0m), is an entangled state in the Unruh quantization: 

|0m) = Y,z~* E{K) \K)l \K)r = J2^ m) S(J,K) \J) l \K) r . (11.28) 

K J,K 

Given that \0m) is an entangled state, when restricted to the R wedge, it fails to be a pure state [14] 
and its description is in terms of a density matrix. The norm of the Minkowski vacuum is now 
given by 

oo ^ 

<M0m> = [I i _ e -2^ m • ( IL29 ) 

m=0 1 C 

The fact that this well behaved state in the plane wave quantization has an infinite norm in the 
Unruh quantization is a consequence of the non unitarily equivalence between these quantizations. 
Recall that different, even not unitarily equivalent, representations of the field algebra are "phys- 
ically equivalent" in the sense of Fell's theorem (see Ref. [5]) and thus one expects to obtain the 
same physical information from the computation of expectation values in either the inertial or 
the Unruh scheme, provided that the Unruh quantization is a faithful representation of the field 
algebra 1 . We can say that our results point in this direction since, as we shall see in Sec. V and 
Appendix C, we obtain, in both quantization schemes, the same physical result. 



III. POSING THE PROBLEM: THE MATRIX OF ENTANGLEMENT 

In this work we shall be concerned with the difference between the initial vacuum state and 
the final state which results from the initial state and the interaction of the classical scalar source 
in uniform acceleration. At the same time we are interested with the restriction of these sates to 
the different regions and thus we shall employ the language of density matrices as indicated by the 
posing of the questions raised in the introduction. 

The trajectory of this source is a branch of an hyperbola lying in one Rindler wedge of spacetime 
which we choose it to be R. The final state of such interaction can be obtained perturbatively by 
the application of the respective S-matrix operator to the inertial vacuum state (this operator will 
be constructed in detail in Sect. IV): 

\f) = S\0 M ). (III.1) 

Recall from Sec. II that the L w modes of the field are zero when evaluated in wedge R and thus 
the scalar source can only excite R modes of the field. Then, for this case the S operator takes the 
form 

S = h ® Sr. (III.2) 



1 We thank Hanno Sahlmann for pointing this out. 
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On the other hand, expressed in the Unruh quantization scheme, state |/} takes the form 

\f)=Y,F(J,K)\J) L \K) R . (111.3) 

J,K 

All the information of the state |/) is encoded in the function F(J,K), and hence one should be 
able to read from the change in this function, AF = F(J, K) — F vac (J, K), the change in the energy 
momentum tensor due to the effects of the source. In particular, this change should, at the same 
time, codify the fact that no change is produced inside the wedge R and that a dramatic change 
occurs in wedge F. Nevertheless, this information is more clearly encoded in the change in the 
density matrix of the state, Ap = Pf — pvac since it can be splitted directly into left, right and 
entangled contributions as we will explain in the following. This is a second reason for using the 
language of density matrices. 

We already know that, when considered as a state of the composite system <g) the 
inertial vacuum state is an entangled state (see [14] for an intuitive explantion of this), and thus 
one expects that also state \ f) will be entangled. In this sense, since both states are pure, their 
density operators cannot be written in the form [15] 

p' = p' L ®p' R , (III.4) 

where p' L and p' R represent the respective partial density matrices defined by 

P L ' R = Tr LiR p. (111.5) 

Nevertheless, one can introduce a traceless operator p e which encodes all the information of the 
entanglement of the state. We propose that the density operators for these states can be written 

as 

p = p L ®p R + (f. (III.6) 

In the Unruh quantization, the p e operator, which we will call the matrix of entanglement, encodes 
the information of the correlation between left and right components of the state. When computing 
expectations of operators localized in either wedges L or R the matrix of entanglement plays no 
role, in fact, it can be shown that for any operator of the form Al®1r and a state described by 
Eq. (III. 6) one has that 

Tr(i L <g> i fl p) = Ti(A L p L ), (III. 7) 

and thus Tv(Al ® lnp e ) = (and the analogous if the operator has the form I/, (g) Ar). The 
information encoded in p e can only be retrieved when computing expectations of operators with 
L and R components and thus, for the case of observables (made up of field operators) this infor- 
mation is only present when one evaluates the expectation values in wedge F. For example, in the 
computations of [4] the information encoded in p e does not enter in the results, which in principle 
are incompatible with the change in {T^ v ) in wedge F. It is therefore that we expected that the 
information about this change would be encoded in p e . 

We can write the inertial vacuum \0m) in the fashion of Eq. (III. 6): 

Pvac = Pvac ® Pvac + Pvac ■ ( IIL 8) 

From Eq. (11.28) it follows that 

Pvac = |0m) (0m\ = Z^EjEk \K) L \K)r L (J\ r(J\ , (III.9) 
JK 
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where we have defined 



Ej = e ~ wE{J) (III. 10) 

(see Eq. (11.24)), and we have introduced the normalization factor Z defined by 

z- 1 ^Y, e2 j ( IIL11 ) 

J 

in order to have Tr(/S„ ac ) = 1. Taking the partial R and L traces in Eq. (III. 9) we obtain respectively 

pt c = ZY J E 2 J \J)ll(J\ , f>tc = Z^ E J \ J )rh(J\ • (IH.12) 
J J 

Using Eqs. (III. 12) and (III. 9) we can write 

Pvac = Pvac ~ Pvac ® Pvac 

= ZJ2EjEk(\J)l \J)rl(K\ r(K\ - ZEjE K \J) L \K) RL (J\ ^K\ ). ( IIL13 ) 

J,K 

For the density matrix for state |/) we have 

PS = Pf ® Pf + Pf- (HI.14) 

Now we turn back to our own specific case and concentrate particularly in the change in the density 
matrices induced by the interaction of the field with the source. Exploiting the fact that for the 
accelerated source the S matrix is an R operator one finds directly from Eq. (III. 1) and Eq. (III. 2) 
that 

Pf = Pvac ® S R ptc S R + Sp £ vac SK (111.15) 

From this last equation we can identify 

Pf=Pvac pf = S R p* ac S R (III. 16) 

and 

p) = Spl ac Sl (111.17) 

We want now to express the change in the state of the field induced by the interaction in terms of 
the change of the density matrix. Note that when restricted to L this change is 5pj = Pf—p^ac = 0- 
The changes 5pJ = pj — p^ ac and Sp e j = p e j — p% ac deserve special attention since they depend on 
the operator S, note that they are traceless. 

The total change in the density matrix can be written as 

5p = pt c ® Spf + bp) % pt c + Spj ® Spf + 5p e 

= pLc®spf + sp e , 

where 5p e = p e ^ — p% ac and we have used dp) = 0. Eq. (III. 18) thus reflects in a clear language what 
we know about the change in the quantum field. As we mentioned in Sect. I, there is theoretical 
evidence that when one restricts state |/) to wedge R one would obtain the same thermal bath as 
that of the inertial vacuum [4], that is, in this wedge there is no effective change in the state. A 
priori one may think that this means that 5pJ = 0. However, the state should have changed in 
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order to account for the radiation emitted by the source, which should be measurable in F. Thus 
one might conclude that all the information of the physical change in the state is encoded in the 
change in the matrix of entanglement, 5p e , which can only be retrieved in wedge F. To be more 
specific, consider an operator of the form A = I/, <g) Ar on &l ® &r. From Eq. (III. 18) one has 
that the change in its expectation is given by 

5Tr(Ap) = Tr(Ap f ) - TV(ip mc ) = Tt(A R Spf). (III. 19) 

Note that the operator Ar can be localized anywhere in Minkowski spacetime as, for example, the 
operator 1 (g) 4>r(x), which is not identically zero outside wedge L. However, if the operator A has 
non trivial L and R components then the change is given by 

STr(Ap) = Tr(i p L vac ® Spf) + Tr(i 5p e ), (111.20) 

that is, in this case the change in the expectation may come from both Spf and Sp e contributions. 

In order to understand the physical change in the state we shall evaluate the change in the 
expectation of the energy-momentum tensor operator T^ u . Recall that in Minkowski spacetime 
(7jjj,(x)) is defined in the point-splitting description as a coincidence limit [5]: 

(f^(x)) f = lunt^,F(x,x'), (111.21) 

X — *x 

where 

F(x, x>) = (/| 4>(x)4>(x') |/) - (0 M | 4>{x)4>(x') |0 m ) (HI.22) 

and is the differential operator 

V = " h^(V a V a ' + m 2 ). (111.23) 

Thus, all one needs to compute the expectation of T^ u is the change in the two point function 
{4>{x)4>{x')) . We have from Eq. (11.20) that 

fa)$(x') = $ L (x)$ L (x') + h{x)^R(x) + $r(x)$ l (x') + Mx)j> R (x'). (111.24) 

To simplify the notation, let us define 

4> = 4>{x), <p' = c/>(x'). (111.25) 

From Eq. (III. 24) and Eq. (III. 18) we have that the total change in the expectation of the two 
point operator is 

Tr(#'«5p) = Tr(^' L p L vac ® 6p R ) + Tr(<^ p L vac «8> Sp R ) + Tr(^' R Sp e ) + Tr(^0' L 5p e ) 

= Tr L (4>J' L rt ac )Tr R (5p R ) + Tr R (4> R 4>' R Sp R ) + Tr{h^p e ) + Tr(4> R 4>> L 5p e ). 

(111.26) 

To get Eq. (111.26) we have used that Tr(i Lji? p e ) = (see Eq. (HI.7)) and that 

Tr L (k pLc) = E E J^ J \ fa \ J ) L = °> ( IIL27 ) 
J 

since states with different number of particles are orthogonal. 

Eq. (III. 26) is the farthest that we can get to reduce Tr((jxj)'5p) using only the fact that we 
are considering a source with support totally contained in R and the properties of the matrix of 
entanglement. Note that the last two terms in the r.h.s. of Eq. (III. 26) are zero when evaluating 
both x, x' 6 L or R since in these wedges the modes and R^ cannot be different from zero 
simultaneously. To go further in our calculation we shall introduce the explicit form of the operator 
S, which we do in the next section. 
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IV. AN ACCELERATED SCALAR SOURCE AND ITS INTERACTION 



We are going to use a scalar source j{x) to model a scalar particle with uniform acceleration. 
Let this scalar current j(x) interact with the field with an interaction Hamiltonian density given 
by 

Hi(x) = ^j(x)4>(x), (IV.l) 

where g is the determinant of the metric. The final state |/) of the field after this interaction is 
given by the application of the S-matrix to the inertial vaccuum state: 



|/) = S|0 M > 5 = fexp 



r^out 

i / H I {x)d 2 x 



(IV.2) 



where T, out are Cauchy hypersurfaces where the interaction begins and ends respectively, and 
T is the time order operator. 

The interaction occurs inside wedge R and thus the in and out Cauchy hypersurfaces should 
bound this region and, at the same time, in order to define states in the Unruh quantization scheme, 
they have to be Cauchy hypersurfaces of the double wedge L U R. We define Sj n as the surface 
constructed by the union of {t = 0, z < 0} and a spatial surface inside wedge R which begins in 
the bifurcation point of the horizons and deviates slightly from the ( = 0, r = — oo horizon (see 
Fig. I). Y, out is defined analogously but its restriction to R is an spatial surface which deviates 
slightly from the £ = 0, r = oo horizon. The initial vacuum state is defined over and the 
final state of the field |/) is defined over T, out . Therefore, although the interaction is present inside 
wedge R, the state |/) is defined after the interaction and one can evaluate expectations in this 
state of operators localized in wedges R and F. 

Our calculation of |/) (and its density matrix) will be in terms of Rindler coordinates for which 
T orders up operators with respect to the time coordinate r. Using Eq. (IV.l) one can put the 
final state of the interaction in the form 

|/)=T(l-i J d 4 Xv^i(x)(A(x)-i J J d^xd^x' ^g^j(x')j(x')^x)4>(x')) \0 M ) + 6(q 3 ), 

(IV.3) 

where the integrations are over the same region as in Eq. (IV.2). It is useful to define the (formal) 
operator 



h = J d 2 x^j(x)<P(x). (IV.4) 

Note that (f>i is of order q. We are going to apply the Wick theorem to the r.h.s. of Eq. (IV.3): 

T(^x)Hx')) = N(fa)4>(xf)) + <0| T^(x)^(x')) |0> , (IV.5) 

where the normal ordering, N, and the vacuum, |0), are with respect to the quantization scheme 
we are dealing with. Using Eq. (IV.5) we expand Eq. (IV.3) to get 

\f) = s |o M ) = (i - sf ) |o M ) - %fa |o M ) - \nhh) |0m) + 0(q 3 ), (IV.6) 



where 



1 roo r-oo ^ 

® = 2 /_ Sx I d2y j{x)j{y) (0mI T \^ x )^y)) 1°^) ( IV - 7 ) 
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is of second order in q. 

A uniformly accelerating scalar particle follows a trajectory which in Rindler coordinates cor- 
responds to the locus of C = Co- From the point of view of Rindler observers, this scalar source 
corresponds to a static scalar current: 

j'(x)=qS(C- Co)- (IV.8) 
This source transforms as a scalar, using Eq. (II. 1) it can be transformed to inertial coordinates 



, w . 9 «*-f+^, (IV..) 

where Co = V a an d « is the proper acceleration of the source. In order to avoid expressions of 
the form x oo one has to introduce some regularization factor to Eq. (IV.8). Inspired in Higuchi 
et.al. [4], who regularize an electric charge and show the consistency of their regulator (we have 
explained their results in Sec. I), we introduce an oscillating factor to Eq. (IV.8): 

j(x)=qcos(er)S(C-Co) (IV.10) 

and at the end of our calculations we shall take the limit 6 — > 0. In fact, for slow oscillations 
{6 <C a) the source is expected to interact with the field as if it were a constant charge q at each r. 

The current j{x) has support totally contained in R, and thus one can see from Eqs. (11.20) 
and (IV.4) that only i?-modes of the field will get excited by the accelerated source. The operator 
4>i takes the form 

h = J R d 2 x^j(x)Mx) = J fa (r u r u + (iv.ii) 

where 

T u = f d 2 x ^j(x)RUx) = <M(w - 6), ^ e = q^smh(ir6)(oK ie (m(o). (IV.12) 
JR 

To get Eq. (IV.12) we have used Eq. (A. 11), lo > 0, and we have chosen to work with 8 > (note 
that the regulator is an even function of 0). Note that Ki ul {z) is real for real uj and z. Finally, Eq. 
(IV.ll) takes the form 

h = g^skihp) Co^(mCo) [h + rj] . (IV.13) 

The role of the regulator we have chosen is to couple the source to an R mode of the field with 
frequency 9 instead of the mode with frequency zero which is somehow pathological. 



V. THE CHANGE OF (f^) AT WEDGES L AND R 

Now we are in a position to compute explicitly 5p R and 5p e in order to evaluate Eq. (III. 26). 
We will work out this calculation perturbatively up to second order in q, which turns to be the 
first relevant contribution. From Eq. (IV.6) we have 

S = (1 - cf) - ifa + Itf^j) + 0(q 3 ), (VI) 
and from this equation we can write the density matrix for the final state as 

Pf = pf+pf+pf+0{q% (V.2) 
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where p^ corresponds to the (second order) renormalized inertial vacuum density matrix given by 

pf = QPvac Q = (l~ 2Re(Sf)). (V.3) 

We shall compute the change in the density matrix with respect to the renormalized vacuum density 
operator: 

8p ren = p f -pf. (V-4) 

Tracing out the left degrees of freedom in Eq. (V.2) we have 

pf = pf {0) +pf {1) +pf i2) +O( q % (V.5) 

where 

pf {1) = -i[h,Pvacl (V.6) 

f>f 2) = \{PvacN{hh) ] +2$ I pKj I + N(]> I 4>l)p*ac)- (V.7) 

H (1) 

Note that pj is traceless. For the entanglement matrix we have 

8? = pf ) +pf ) +0{q\ (V.8) 

where and p e ^ are defined as Eqs. (V.6) and (V.7) changing R — > e. 
The first order contribution to Eq. (III. 26) is 

Tr(M>'SpW) = Tr L (t L fi L rt ac )Tr(pfW) + Tr(<^p* (1) ) + Tr{(M R + &fo)pj (1) ). (V.9) 

The first term in the r.h.s. of Eq. (V.9) is zero since Pf is traceless. From the expression for 
Pf i Eq. (V.6), it can be proved that 

Tt{4> R $ R pf) = Tr{4>J R pf) = Tr(M' L pf ] ) = 0. (V.10) 

This can be understood heuristically from the fact that these traces represent a sum of brackets in 
the Minkowski vacuum of three field operators, which are necessarily null. From Eq. (V.10) and 
the definition of p^ we have 

Tr(<^p; (1) ) = T±{4> L $Rpf) ~ ^{h^RPvac ® pf 

= M4>l4>' r p?) (V- 11 ) 

= 0. 

To get the second equality in Eq. (V.ll) we have used Eq. (III. 27) to conclude that 

T±(h$R Pvac ® pf W ) = Tr L (4>lpLc)^r (4>' R pf {1) ) = . (V.12) 

Note that this last equation is valid for any order. We have then proved that Tr(<f>(f>'8p( 1 ')) = 0. 
Now we turn to the second order contribution, 

Tr(#'<5p( 2 )) = Tr L (<^' L p£ ac )Tr(pf 2) ) + Trfafof™) + Tr((M' R + M R )pf ] ) ■ (V.13) 
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For the mixed LR terms, we have that, similarly to Eq. (V.ll) 

Tr(j> L $' R pf ] ) = T±{h$R pf) ~ T±{M' R Pvac ® pf (2) ; 
= Tr(4>J R pf). 



(V.14) 



Using Eq. (V.7) it can be proven directly that (see Appendix B) 

T^L&pJ (2) ) = ZV 2 e Tr(4> L 4>' R p vac ), (V.15) 
Tr(M' L p e / 2) ) = ZV 2 e Tr(4> R 4>' L p vac ) . (V.16) 

Also from Eq. (V.7) (see its basis expansion in Eq. (B.10)) we have that Tr(p^ 2 ^) = ty 2 , and thus 

TrUhftLPvac)Tr(pf (2) ) = Z^ 2 Tt l {^lpLc) • (V.17) 

These terms, since are proportional to the corresponding LR and LL parts of the two point function 
in the vacuum, are going to be absorbed in the renormalization of the final result. We have then 
obtained the, in principle, unexpected result that at least for the change in the expectation of the 
Tpj, operator there is no contribution from the change in the entanglement matrix, although as one 
can see from Eq. (III. 17) there is actually a change in this operator. All the information of this 
change comes from the change in the R density matrix of the state, SpJ. As shown explicitly in 
Appendix B we have that 

Tr(4> R 4>' R pf 2) ) = W 2 lm[R e (x)} lm[R 6 (x')} + Z^ 2 Ti(p^ c 4> R (x)4> R (x')) . (V.18) 

Finally, adding up Eqs. (V.15)-(V.18) we have that the change in the two point function between 
the inertial vacuum and the state generated by the interaction of the scalar source is 

Tv(4>(x)4>(x')6p ren ) = W 2 lm[Re(x)} lm[R e (x')} + ^Tx{^x)4>{x')pf) + 0(q 3 ), (V.19) 

which is valid for all x, x' £ M and we have used Eq. (V.3) to express p vac in terms of the 
renormalized vacuum density operator p^p. However, as we have said, we still have to renormalize 
this change in the expectation value. We define the renormalized field operator as 

*2xl/2 



<f> ren (x) = I 1 - I <f>(x), (V.20) 

and thus the renormalized change in the expectation of the two point function is 

C e (x, x') = Tr(4> ren (x)4> ren (x')6p ren ) = W 2 lm[R e (x)} Im[^(x')] + 0(q 3 ) . (V.21) 

The second order term in this change reads 

Cf\x,x') = 4V 2 9 Im[Re(x)] lm[R e (x')} . (V.22) 

We have arrived to a regular expression for the change in the expectation value of the two point 
function and thus now we can reconsider the source as static by taking out the regulator taking 
the limit 9^0. In the following we are going to evaluate explicitly Eq. (V.22) in wedges R and F. 
First, note that form Eq. (IV. 12) and the fact that Bessel functions Kig(mCo) are regular whenever 
Co / we have that 

lim ^ e = . (V.23) 

0->O 
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Note that ^ is an overall factor in cf\x,x'). When computing the expectation of 4>{x)(f)(x') the 
correct procedure is to take the limit 6 — > at the end of the calculation since the field operators 
couple to the frequency 6 of particles in |/) which are emitted and absorbed by the source. 

First we will work the case when x, x' G R, for we can express these points in terms of Rindler 
coordinates: x = (r, Q, x' = (t', £'). Unruh modes R^, when restricted to wedge R take the form 
of Eq. (II. 4) (see also Eq. (A. 11) in the Appendix). Using the fact that functions Ki u) {z) are real 
whenever to is real and z > 0, we have that 

lm(R e (x)) = _ V smh (^) sin (fl r )jQ e ( m ^ x£R _ (v.24) 

7T 

And thus, from Eqs. (V.22) and (IV. 12) 

Cf\x,x') = ^q 2 smh 2 (Tr9)($K? 6 (m(o) sm(8T) sm(8T')K ie (m()K ie (mC) x,x' G R. (V.25) 

Recall that R is an open set bounded by the horizons and thus £, ^ 0, so the Bessel functions 
Kig(mQ, Kig(m(') are regular. Then we have that the second order change in the two point 
function between the inertial vacuum and state |/) in wedge R is 

lim Cf ] (x, x') = x, x' G R . (V.26) 

This resut, is consistent with the fact that, as pointed out in Ref. [4] the source is in thermal 
equilibrium with the field inside wedge R. Any observer inside this wedge will not be able to notice 
any change in the expected value of T^ v due to the presence of the accelerating source. Note that 
the specific form of the two point operator has played a crucial role to get to Eq. (V.22). 

Now we proceed to evaluate Eq. (V.22) in wedge F. Recall that in F, r is a spatial coordinate 
and C is timelike (cf. Eq. (A.2b)). Form Eq. (A. 12) we have that the mode Ro{x) restricted to F 
takes the form 

R oW = ~W2 f. Z m [^SVO + e-^SVC)] x e F, (V.27) 
2°' z A /2smh(7r6 l ) L J 

where coordinates (r, £) are defined in Eq. (A.2b) and H^ ),{2) are the first and second Hankel 
functions. Using the definitions of H$'^ in terms of Bessel functions, Eqs. (A. 19) and (A. 20), it 
follows immediately that 

lm(R 9 (x))=- 1 = (e- ieT J- i0 (mO + e idT J if) (mO) x G F . (V.28) 
4y / smh(7r#) 

Now using Eqs. (IV. 12) and (V.28) we have that 

lim Cf\x, x') = q 2 (%Ki(m(o)Mm()J (m(') for all x,x' G F . (V.29) 
o — >o 

Compare with Eq. (V.26). This expression is the (non zero) change in the two point function in 
wedge F due to the interaction. It contains the information of the field radiated away from the 
source into F. It should be noted that Eq. (V.29) is not valid at the horizons. In Appendix C we 
make the computation, in an inertial frame, of the same change in the expectation of the two point 
function, Tr((j)(x)(j)(x')5p), for x,x' G F and obtain exactly Eq. (V.29). Therefore, at least for the 
particular case we are dealing with, both quantum descriptions produce the same physical results. 
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VI. AN EXAMPLE OF LEFT-RIGHT INTERFERENCE: TWO SOURCES WITH OPPO- 
SITE ACCELERATION 

As we have seen in the last section, the matrix of entanglement plays no significant role in 
the change in the expectation of T^ v for the interaction of an accelerated source with the inertial 
vacuum. However, the presence of another (accelerating) source lying in L may affect the entan- 
glement of the final state and thus p e may play a part in the change of the expectation of T^ u . In 
this section we work out these calculations. 

Suppose that additionally to the scalar particle considered in the latter sections we have an 
extra scalar particle in uniform acceleration in wedge L with scalar charge q. The perturbation of 
the field due to this pair of particles is given by 

(vu) 

where 

j L (x) = qcos(eT L )5(( L - Co) j R (x) = qcos(eT R )5(( R - Co) • (VI.2) 

Rindler coordinates in wedge L, (tl,(l) are given by Eq. (A. 2a). To avoid confusion, all along 
this section we will call Rindler coordinates in R: (t r ,(r) (originally we have used (r, £)). Note 
that in order to have independent sources we have introduced the cosine regulator with a different 
parameter for the source in wedge L. 

Using Eq. (IV. 1) for j{x) given by Eq. (VI. 1) we have that 

S' = fexp [ - i[(4>f ® + (1 L ® 4>?)}] = f exp[-itf] ® exp[-i$f] = S L ® S R , (VI.3) 

where 

ft = * § {l~ e + /]) , $ = *e {re + r\) , (VI.4) 

and T is the time order operator. In this case we have two different time parameters, tx and tr 
and thus T will time order L and R operators independently (recall that 4>l and 4> R commute). 
The factor is given by Eq. (IV. 12) and ^ is now 

m~ e = q^smh(n9)( K i§ (m( ) . (VI.5) 

Let \g) be the final state of this interaction, analogously to Eq. (III. 15) now we have that its 
density matrix takes the form 

P 9 =Pg®f>f + Pg, (VL6) 

where 

p L g = S L p L vac Sl pf = S R p* ac S R , p e g = S'pl ac S'K (VI.7) 
Then we have that the change in the density operators is given by 

5pg = Pvac ® Spf + 8p L g ® p* c + 5p L g ® 5pf + 5p e g , (VI.8) 
where 5p g = p g — p vac and all other differences are defined analogously. 
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Now we compute the change Tr (<fi(j)' 5 p g ), analogously to the case of the single accelerating 
particle the first relevant term is of second order in q and q. Using Eq. (III. 24) and Eq. (VI. 8) It 
can be seen that 

Tr(<^/f )) =Tr L (fa&pj;™) + Tr R (M' R 6pfW) 
+ Tr((4> L $ R + 4> R 4>' L ) 5p L g V ® ^W) 

The first four terms in the r.h.s. of Eq. (VI. 9) are analogous to Eq. (V.17) and and Eq. (V.18) for 
the single source case. In effect, we have that 

Tr(^tof 2 )) = ^ 2 e Im[Re(x)] lm[R e (x')] + Z^Tr(p^ c 4> R (x)Mx')) , (VI.10) 

Tr(4> L 4>' L 5pW) = Wpm[L § (x)] lm[L § (x')] + ZV}Tr(tf a J L (x)j> L (x')) . (VI.ll) 
Similarly to Eq. (V.6), it can be seen that 

8p L g V = -i^ptc-Pvactt), (VL12) 

8pfW = -i^fp^-pKjf), (VI.13) 

and from these equations we have that 

Tr((4> L $' R + 4> R 4>' L ) 5p L g ^ ® SffW) = W § * (lm(L § (x))Im(R (x>)) + lm(R § (x))Im(L (x'))) . 

(VI.14) 

On the other hand, it can be shown that the second order contribution from the matrix of entan- 
glement is 

Tt((4> l 4>' r + 4>R4>' L )5p g W) = Z(*j + *§) Tr((4> L 4>' R + fafopvac), (VI.15) 

which, as in the case of a single source, also corresponds to terms that will be absorbed in the 
renormalization of the final change in the expectation value. Adding up all the contributions, 
Eq. (VI. 9) takes the form 

Tr^/f )) =4[* § 1m(L § (x)) + *ehn(Ro(x))] [* § 1m(L § (x')) + ^gIm(R e (x'))] 
+ Z(y} + yj)Tr$(x)4>(x')p vac ). 

Thus, the renormalized change in the two point function reads, up to second order 

Tr(0 re J' ren 5p gren ) = 4[^ § lm(L § (x))+^ e lm(R e (x))] [* § Tm(L § (x'))+*eTm(Rg(x'))]+. . . (VI.17) 

After taking the limits 6,9^0 and evaluating at x = (tf, Cf) and x' = (t' f , £' f ) we obtain 

Jim Tr((f> ren (x)(i> ren (x')5pg ren ) = (q(oK (rn( ) + q(oK (m(o)) 2 J (m( F )J (m( F ) . (VI. 18) 

0,0-»O 

As expected, if we turn off the charge in wedge L (q — ► 0) we recover our previos result, Eq. (V.29). 
For the case of the two accelerating sources it turns out that also all the contribution to the change 
in the expectation value of the two point function comes solely from the change in the partial 
matrices, Sp g , 5f>^- In particular, the interference term, Eq. (VI.14), is determined by the latter 
pair of matrices, that is, for the case we have just analyzed all the information of the change in 
the expectation of T^ v is only encoded in 5p g , 5p^. 
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VII. DISCUSSION 



One of the main goals of our work was to reconcile the fact that the final state of the field 
appears to remain an undisrupted thermal state in both the left and right Rindler wedges, with 
the expected change induced by the source on field observables, such as the energy momentum, in 
the future wedge. This issue is recasted in terms of the codification in the state of the field of the 
pertinent information that exhibits the change in the expectation value of T^ v . 

At the beginning it was our belief that, since there is no change in the expectation of T^ v 
in wedge R (neither in L) the physical change in wedge F cannot be induced by the particular 
behaviours of the state when restricted to either wedge. Hence, the information of this change 
information should have been encoded in some part of the state which is not represented by any of 
the restricted density operators pj and /5j (or in their respective changes). It is in this sense is that 
we proposed the decomposition of pf given by Eq. (III. 6) with the particular introduction of the 
matrix of entanglement p e . As we explained in Sec. Ill this operator plays no role when computing 
expectations of observables localized in wedges L and R and thus was a good candidate to account 
for the change of {T^ v ) in wedge F. Nevertheless, we computed this change perturbatively and 
found that it has contributions only from the change in the density matrix describing the state 
in the wedge R (see Eqs. (V.18) and (V.21)). That is, when evaluating (T^ u ) in either wedges R 
and F, its change is determined solely by the characterization of the state in wedge R. This result 
contrasts with our initial expectation that the information about the change of {T^ v ) would be 
encoded in the change in the entanglement matrix, 5p e . 

In order to obtain this result, we had to introduce by hand a particular regularizing function, 
cos(0r), into the current describing a uniformly accelerating scalar source with the prescription to 
take the limit — ► at the end of the calculation. As mentioned earlier, this regulator was inspired 
by a similar computation done in [4]. Despite the seemingly artificial choice of this regulator, we 
have shown that when making the regulator independent calculation of the change in the same 
expectation {cf)(x)4>(x')) , using a plane wave quantization scheme, one obtains the same results as in 
the Unruh scheme with such regulator (see Appendix C). One can give an heuristic explanation to 
the fact that this regulator is physically correct as follows. In principle, if one is precise, one would 
like to describe the radiation due to a real physical particle, which should be described as a quantum 
object itself. However, this description has a serious drawback regarding our wish to describe 
the source as a uniformly accelerating particle: A quantum particle does not move in a definite 
trajectory and thus assigning to it a particular acceleration is impossible. On the other hand, the 
nature of the quantum field is distributional and therefore the correct quantum description of the 
interacting source should be in terms of test functions of compact support (note that Eq. (IV. 10) 
lacks this property) and therefore, the source would correspond to an extended object. To such 
object we can not naturally ascribe a uniform acceleration: if the object is to maintain "its shape" 
along its trajectory then different parts should have different proper accelerations. However, we 
know from [4] that a treatment using classical point-like sources (with definite proper acceleration) 
together with a certain type of regulator produce physically correct results (in particular, results 
that are fully consistent with the Equivalence Principle). 

The regulator used in [4] consisted in the introduction of an artificial oscillation with frequency 
9 in the strength of the source, to identify therefore expressions of the form Ox oo occurring in the 
calculation and to proceed to carry all calculations to the end before taking the limit 9 — > 0. We 
are thus assuming that the introduction of such regulator, along with the prescription to take the 
limit 9 — ► at the end of the calculations reflects in an effective way the description of a quantum 
source in uniform acceleration interacting with the field. (Nevertheless, the robustness of the result 
would be ensured if one confirms that the same physical behaviour is to be obtained for a wider 
class of regularizing functions.) It is however worthy to emphasize that in the inertial calculation 
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of Appendix C there was no need to introduce such a regulator. 

From these considerations it follows that our calculations make sense only if the limit 8 — ► 
is taken at the end of the computation of the expectation values. Actually, one may be tempted 
to take this limit directly into the density matrix pf of the state. Ignoring for a moment details 
regarding the precise notion of limit of an operator, one can see that due to the overall factor vl/g, 
which goes to zero as 9 — > 0, every single term in the expansion of pf is also zero except for those 
terms proportional to p vac - Thus, one would conclude that there has been no change in the state 
of the field in R due to the presence of the source. As this was the only potential contribution to 
the change in the expectation of T^ u , we would be led to the erroneous conclusion that there is no 
change in this quantity. 

From the calculations in Sect. V we concluded that the information of the change in (T^ u ) is 
encoded in Sp R . We now want to consider how is it encoded. The answer to this question relies 
on a subtle interplay of the field operator (p(x) and the density matrix in the present formalism. 
Let us focus on the details of our specific calculation: The overall factor ^ in the second order 
contribution to pj comes from the fact that the source is located in wedge R. In this wedge Unruh 

modes Rg(x) — > as 6 — > (whenever x is not at the horizon). When computing Tr(4>R4> R pjr ) 
the field operators are sensitive to the frequency 6. In fact, they only excite modes with uj = 9 as 
could be expected on Rindler energy conservation grounds. The particular form of pFj" determines 
the structure of the contribution given by Eq. (V.22) which in turn, due to the different behaviors 
of the Unruh modes in wedges R and F, is zero in the former and not zero in the latter. 

We interpret these results by saying that the "0" Rindler energy modes, which were of such 
concern in regard to the definition of the theory (see discussion at the end of Sec. II A), are in 
effect, essential in order to obtain in the accelerated frame description identical results as in the 
inertial one. Physically we could think that these modes are excited by the slightest quantum 
fluctuations of a realistic quantum particle and that their excitation would be directly felt in the 
future wedge. These results seem to be in accordance with the spirit of those obtained in [16] 
where it is argued that the zero energy modes seemed to be undetectable (with an appropriate 
definition of detectability) when confining the detection to the right wedge. Finally, it is our belief 
that this work has helped in clarifying the questions raised at the beginning. Furthermore, there 
is a technically analogous situation which indicates that a stationary particle just outside of the 
horizon of a stationary black hole could be "emitting" towards its interior. This work shows a clear 
path to studying the changes in the energy momentum tensor in that situation. 
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Appendix A: UNRUH MODES AND REPRESENTATIONS OF BOOST MODES 

One can generalize Rindler coordinates to all Rindler wedges assigning always the respective 
coordinate r to the parameter associated to the generator of boosts about the origin in the z 
direction: 




(A.l) 
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In wedge L one should do an exception in order to to guarantee that the future direction coincides 
with that of inertial time, in this case, r M is — b^. In wedges F and P the coordinate r is spacelike. 
We have [17]: 

z = — Ccosh(ar) i = £sinh(ar) (t, z) G L, (A. 2a) 

z = Csinh(ar) t = ( cosh(ar) (t, z) G F, (A.2b) 

z = - C sinh(ar) t = -( cosh(ar) (t, z) G P, (A. 2c) 

where £ > in each region. 

Boost modes, Eq. (11.13), should be thought of as distributions and thus one cannot evaluate 
them in one particular point. However, what we can do is to apply this distributions to test 
functions which have support defined on a certain open region. The integral which defines the 
boost modes can be expressed in terms of the accelerated coordinates (r, Q given by Eqs. (II. 1) 
and (A. 2). As Unruh modes are defined by boost modes, from this operation we can express Unruh 
modes in accelerated coordinates too. 

For example, let / be a test function in M with supp(/) C R, the evaluation of B u at function 
/ is given by 

Bjf] = — 1- r d9 e-' lUj9 f d 2 X e im(zsmh(e)-tcosh(e)) f ^ z , _ ( A 3 ) 

2 d / 2 7r y_oo J R 

Since the space-time integration is in wedge R we can express the d 2 x integral in terms of Rindler 
coordinates (c.f. Eq.(II.l)): 

B ^ = 23?^ d0 e ~ W(> J R d2% e tmCsinHe - T) f(r, C) . (A.4) 
Changing the integration order and using the following relation for the Bessel functions [18] 

K v (x) = -e*™ / da e ^ e -^ sinh ( a ), (A.5) 

2 J~oo 

it can be seen that 

B u [f] = -±j=e"f f d 2 xe-^K iui (mOf(r,0 • (A.6) 

7TV2 JR 

By analogous arguments it can be seen that in other wedges boost modes take the following form 
[13] (coordinates in the following eqs. are respective to each wedge defined by Eqs. (II. 1) and 
(A.2)). 

BMt,Q = -^e^e-^K^mC), (A.7) 
B u \l(t,Q = -^e^e^KiUmQ, (A.8) 

BMt,C) = -^efe-^K), (A.9) 
BMr,C) = ^e^e-^K), (A.10) 

where H^'^ are Hankel functions. Note that in wedge L, B^{t, Q has negative frequency w.r.t. r 
because we have chosen the time translation generator in L to be = —b^, where b^ is the boost 
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generator. Using Eqs. (A.7-A.10) in Eqs. (11.18) one obtais directly the representations of Unruh 
modes in each wedge. Modes R w are given by: 



Ru\l{t,Q 

RMt,C) 



i Vsinhvrwe-^^ (m() , 



2 3 / 2 v /2sinh(vr^) L 



= 0, 



And La, modes are 



Loj\r(,t,Q 
Lw\f(t,() 



2 3 / 2 V2^mh 

o, 



2 3 / 2 ^2sinh(7rcj) 



1 



7T 



C) = -Vsinhvrwe lu}T K iuJ (m(), 
L u \p(t,Q 



2 3 / 2 V2ihrh 



^(mO + ^VO 



(A.ll) 

(A.12) 
(A.13) 
(A.14) 

(A.15) 
(A.16) 

(A.17) 
(A.18) 



Here we put some useful relations for the Hankel functions [18] from which one can simplify the 
expressions for Unruh modes 



Hi 2 \z) = 



1 



i sin(i/7r) 
1 



i sin(^7r) 

These latter equations are used to obtain Eq. (V.28). 



[j^(z)-e-^J u (z)], 
[e^Mz) - J_„(z)] . 



(A.19) 
(A.20) 



Appendix B: SECOND ORDER CALCULATIONS 

Along this work we do several second order calculations of expectation values of <j>{x)<j)(x') . In 
this Appendix we want to put the details of the calculus of Tr(/5j <j)R{x)(j)R{x')) which leads to 
Eq. (V.18); all the other calculi are analogous to this one. 

The state |/) takes the form of Eq. (11.26) in the Unruh quantization, in effect, from Eq. (IV.6) 
one can express |/) as (from now on we will omit the (8> in the Unruh states) 

\f) =J2 F (J,K)\J)l\K) r , F(J,K) = F (J,K) + F!(J,K) + F 2 (J,K) +0(q 3 ), (B.l) 
J,K 

where the first term is Fq(J,K) = QF vac (J, K) (Q is defined above Eq. (IV. 7)), and F vac (J,K) is 
the left-right superposition function defining the inertial vacuum state \0m) (cf- Eq. (11.28)): 



F vac (J,K)=Z 1 / 2 e-* E ^5(J,K), 



(B.2) 



where Z is the normalization factor defined in Eq. (III. 11). In order to express the other terms in 
Eq. (B.l), it is useful to define the normalization factor Nq(K) by 



N?(K) = 



VK e + 1 a = + 
\[Kq a = — 



(B.3) 
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where Kg is the particle component of state \K) L R in the mode with frequency centered on to = 9. 
Using Eq. (IV. 13) in Eq. (IV.6) we find 

F {J,K) = QZ^e'^^S^K), (B.4) 
F ± {J,K) = -izVH B e-* E ( J) (N+{J)5{K, J + lg) + N e (J)5{K, J-l e )), (B.5) 

F 2 (J,K) = ~Z 1 / 2 *le-^(N e (J)N e (J-l e )5(K,J-2e) 

+2N e (J)N+(J - l e ) 6(K, J) + N+(J)N+(J + l e ) 6(K, J + 2 fl )) . (B.6) 

Eq. (B.4) is actually Eq. (11.28) with the extra Q factor; E{J) is defined by Eq. (11.24). The 5 
functions are defined by Eq. (11.25). The J — lg which appears in the last term in the r.h.s. of 
Eq. (B.5) corresponds to a normalized state | J — 1$)r defined by 

rg\J)R = N 9 {J)\J -l e ) R (B.7) 

with particle content J — lg = {J Uo , . . . , Jg — 1, . . . }. The other terms are defined analogously. 
The density matrix of state |/) reads 

Pf = \f)(f\ = E F(J,K)F(j',K'y \J) L \K) RL (j'\x(K'\. (B.8) 

J,K,J',K' 

Let's write it in the following manner 

p f = pf + pf + pf+0^). (B.9) 
Taking the trace over the left degrees of freedom to this expression we obtain that 



pf 2) = Z^g 



\{e-™-l?p a -\{e™-l?p b + Pc (B.10) 



where 



Pa 


-E« 




K 


Pb 


= £e 




K 


Pc 


-E« 

K 



(B.ll) 
(B.12) 
(B.13) 



Writing the field operator as 



te= E E^wc (B.i4) 

ct=+,— m=0 

where 

#^ m ( x ) = K m ( x ), R Z m ( x ) = Ru> m {x), r+ m = rl m , f~ m = f Wm , (B.15) 

we have 

oo 

Tr(pJ R (x)4> R (x>)) =J2Y, E ^»V^v^TI<W^(if + ll)x 

a, a' m,n=0 

x KJx)K'Jx')r(k + 2 e \K + lZ m + , (B.16) 
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where \K + lj m )i? = \K ± \ Wm ) R . The expectation value appearing in the r.h.s. of Eq. (B.16) gives 



R\K + MK + lZ m + lZ n ) R = <W <*«',+ <W S Unfi ■ (B.17) 

From this one can read that there is only a contribution to Eq. (B.16) when the operator 
4>r{x)4>r{x') creates two particles in the mode u = 6 in the state defined by p a . Using Eq. (B.17) 
in Eq. (B.16) we have that 



Tr(pJ R (x)Mx')) = R* e {x)Rl{x') Y e- 2 * E{K \Ke + 2)(K e + 1) . (B.18) 

K 

To evaluate the sum in the r.h.s. of Eq. (B.18) note that from Eq. (11.27) we have that 

oo oo oo 

Y j e- 2 * E ^f{K e )= Y e- 2 * eKe f(K e )x [] Y e -^ m K^ m 



K 



OO 



(B.19) 



(0m|0m) (l-e^ e ) Y, e- 2 ^f(K e ), 

Kg=0 



where f(Ke) is any function of K$ and (0m|0m) is given by Eq. (11.29). For latter use it is 
convenient to define 

oo 

Ge[f(K 6 )] = Y ^ 2neKe f{K e ) . (B.20) 

Kg=0 



Using Eq. (B.19) then Eq. (B.18) becomes 



Ty( P J r {x)^ r {x')) = (0 M \0 M ) *_ 27t6 2 R* e (x)R* d (x'), (B.21) 



where we have used that 



G e [(K e + 2)(K e + 1)] = _ g 3 . (B.22) 

Analogously we have 

Tr(p b 4> R (x)4> R (x')) = (0 M \0 M ) 2 J_ 2 Re{x)Re{x') . (B.23) 

Before computing Ti( y p c <p(x)(j)(x')) let's define 

H Um [K\(x,a/) = K^KJx^ix') + (K^ + ^R^R^Jx') . (B.24) 
It can be verified that 

oo 

Tr{p^Mx)Mx')) =Y^ E{K) Y H« m [K](x,x') . (B.25) 

K m=0 
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Now we compute 

Tr(^e- 2 ^W [e 2 * e K e + e- 2 * e (K e + l)\ \K) RR (K\4> R (x)4> R (x')) = 

K 

=J2 ^ 2wE(K) [{^ d K 2 e + e-™Ke{K e + l))R* e (x)R e (x') 

+ (e 2 ^K e (K e + 1) + e~ 2 ^(K e + l) 2 )R e (x)R* e (x')] 

oo 

+ J2e- 2 ^ K \e 2 ^K e + e- 2 ^(K e + l)) £ H^[K](x,x>) 

K m=0 



K 



(B.26) 



From the definition of G$[f(Kg)], Eq. (B.20), it can be shown that 

e 2 ^G e [K 2 ] + e- 2 * e Ge[K e (K e + 1)] = G e [K (2K + 1)] + Gg[l], (B.27) 
e 2 ^G e [K e (^ + 1)] + e -^ e G e [(K e + l) 2 ] = G e [(2K e + l)(K e + 1)] + G g [l], (B.28) 
e 2 ^G e + e - 2we G e [(#«, + 1)] = G e [(2K + 1)] . (B.29) 

Now we use Eq. (B.19) to simplify the r.h.s. of Eq. (B.26). From Eqs.(B.27)-(B.29) we have that 
Tr (E e ~ 27r£W [e 27ve Ke + e- 2 ^(K e + l)] \K) RR (K\4> R (x)4> R (x')) = 

K 

oo 

= (0m|0m) {R* e (x)R e (x') + R e (x)R* e (x')) + ^e" 2 ^W(2^ + 1)) ^ HujK](x,x') . (B.30) 

K m=0 

Then, from Eq. (B.30) we have at once that (see Eq. (B.13)) 

Tr(p c Mx)4> R (x')) = (0 M |0 M ) (RU^)Re(x') + Re(x)Rl(x')) +Tr(p^ R (x)4> R (x')), (B.31) 
where we have used Eq. (B.25). Finally, from Eqs. (B.21), (B.23), (B.31) and (B.10) we have 

Tr(pf 2) Mx)Mx')) =ttg( - R d (x)R 6 (x') - R* e {x)R* e {x') + R* e {x)R e {x') + Re{x)R* e {x')) 

+ Zy 2 Tr{pK c 4> R (x)4> R (x'))- 
and Eq. (V.18) follows directly. 

Appendix C: INERTIAL CALCULATION 

In this Appendix we want to show that one obtains exactly Eq. (V.29) when computing the 
change in {(p(x)4>(x')) in the inertial scheme. For this case we are going to use the scalar source 
given by Eq. (IV. 9) written in the form 



j( x ) = iCo 1 a u , (C.l) 

V 1 + Co 



where (o = 1/a and a is the acceleration of the source. Using the inertial field operator 

/oo 
dp (i> P (x)a p + ipp(x)al) , (C.2) 
-oo 
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where ip p (x) represents a plane wave with frequency uo p = +\/p 2 + m 2 , 

M x ) = r~ l ; e-^" t+ipz , (C.3) 

in Eqs. (IV. 1) and (IV. 2), it can be seen that the second order renormalized change in {(p(x)(p(x')) 
between states |/) and \0m) is given by 



CTJ(x,x') = 4Im[Q(x)]Im[Q(x')], (C.4) 

where 

r f'OO 

Q(x)= d 2 x' dpj(x)r p (x')Mx)- (C5) 

J M J — oo 

Recall that for the quantized scalar field we are considering the positive frequency function is given 
by 

/oo 
d P r p (x')M x ) (C6) 
-oo 

and from [6] we have the following result 



AW(x,0) = - x < 



g| 2) (mV ^?) t>\z\ 

%K (mVz 2 - 1 2 ) \t\ < \z\ . (C.7) 
-H^imVt 2 -z 2 ) t<-\z\ 



From now on we shall suppose that x G F. Using the fact that AW(i,i') = AW( x — x', 0) and 
Eq. (C.7) it can be seen that 

Im[Q(x)] = ^ f 7 =^== J (m^MF)), (C.8) 



where 



a(0 = i (-(t - t') 2 + {z- ^' 2 + C 2 ) 2 ) (C.9) 

/2 s ) 

and cr(i _ ) = 0. To get Eq. (C.8) we have used that Hi, (y) = J u {y) — iY v (y) and the fact that Jo(y), 
Yo(y) are real for y > 0. Making the changes of variables t' = Cosinh(r/Co) and u = y / —2a(r) we 
obtain 

lm[Q{x)] = q -^ T uJ ^ mU \ ) ^=du = q -^K (mto)MmO, (CIO) 

2 y v/4C 2 C 2 + (u 2 + Co 2 - C 2 ) 2 2 

were we have used x = (t,z) and i = Ccosh(r p /Co), z = Csinh(r p /Co) (see Eq. (A. 2b)). The 
derivation that leads to the second equality in Eq. (C.10) is analogous to that which leads to Eq. 
§13.54(1) of [11]. Thus we have proved that 

Ci 2) (x,x') = q 2 ( 2 K (mCo) 2 MmOMm('), (C.ll) 

which coincides functionally with Eq. (V.29). In this computation one has to apply Wick theorem 
with the notion of time and normal ordering associated to the inertial time parameter t. However, 
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the effect of this choice of time does not show up in the physical change of the two point function 
but only in renormalization terms. Note that, in contrast to the accelerated frame calculation, we 
did not need to introduce any regulator into the current (neither any cutoff as in the inertial frame 
calculation in [16]). 
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